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ABSTRACT 

We  have  developed  a  dynamic  geomagnetic  vertical  cutoff  rigidity  model  that  predicts  the  energetic  charged 
particle  transmission  through  the  magnetosphere.  Initially  developed  for  space  applications,  we  demonstrate 
the  applicability  of  this  library  of  cutoff  rigidity  models  for  computing  aircraft  radiation  dose.  The  world  grids 
of  vertical  cutoff  rigidities  were  obtained  by  particle  trajectory  tracing  in  a  magnetospheric  model.  This  refer¬ 
ence  set  of  world  grids  of  vertical  cutoff  rigidities  calculated  for  satellite  altitudes  covers  all  magnetic  activity 
levels  from  super  quiet  to  extremely  disturbed  (i.e.,  Kp  indices  ranging  from  0  to  9^)  for  every  three  hours  in 
universal  time.  We  utilize  the  Mcllwain  "L"  parameter  as  the  basis  of  the  interpolation  technique  to  reduce 
these  initial  satellite  altitude  vertical  cutoff  rigidities  to  cutoff  rigidity  values  at  aircraft  altitudes.  Published  by 
Elsevier  Ltd  on  behalf  of  COSPAR. 


INTRODUCTION 

We  have  a  complete  set  of  world  grids  of  vertical  cutoff  rigidities  each  5°  in  latitude  and  5°  in  longitude 
for  a  spacecraft  orbiting  at  450  km.  These  vertical  cutoff  rigidities  were  obtained  by  particle  trajectory  tracing 
in  model  magnetospheres  and  include  values  for  all  magnetic  activity  levels  from  super  quiet  to  extremely  dis¬ 
turbed  (i.e.,  Kp  indices  ranging  from  0  to  9"^)  for  every  three  hours  in  universal  time.  We  used  the  Tsyganenko 
(1989)  magnetospheric  field  model  combined  with  the  International  Geomagnetic  Reference  Field  for  Epoch 
1995.0  (Sabaka  et  al,  1997).  The  Tsyganenko  (1989)  magnetospheric  field  model  describes  the  magneto¬ 
spheric  field  topologies  for  Kp**  magnetic  indices  from  0  to  5.  The  Boberg  et  al.,  (1995)  extension  was  used  to 
describe  the  magnetospheric  fields  for  magnetic  activity  levels  exceeding  Kp  values  of  5.  For  convenience  we 
have  labeled  these  as  Kp  6  through  10  in  Dst^  increments  of  -100  nT. 

CUTOFF  RIGIDITY  DETERMINATION  PROCEDURE 

Cosmic  ray  trajectory  calculations  were  initiated  in  the  vertical  direction  from  a  distance  of  6821.2  km 
from  the  geocentric  (450  km  altitude  above  the  average  earth  radius  of  6371.2  km).  We  utilized  the  Bulirsch- 
Stoer  numerical  integration  technique  (Press  et  al.,  1989);  each  step  length  was  about  1%  of  the  gyro-distance. 
The  "sensible"  atmosphere  of  the  earth  was  considered  to  extend  20  km  above  the  international  reference  ellip- 


*  Rigidity  is  momentum  per  unit  charge  and  is  a  canonical  unit  that  is  especially  useful  in  characterizing  charged  particle 
access  in  magnetic  fields.  All  particles  having  the  same  magnetic  rigidity  will  have  identical  trajectories  in  the  magnetic 
field,  independent  of  elemental  or  isotopic  composition,  particle  mass  or  atomic  charge. 

**  Kp  is  a  quasi-logarithmic  magnetic  activity  index  that  ranges  from  zero  to  9  in  major  increments.  There  are,  3  sub¬ 
increment  levels  (-,  0,  and  +)  within  each  major  increment. 

^  Dst  is  an  index  of  the  magnetic  fields  generated  by  during  magnetic  storms,  measured  in  units  of  Nt  at  the  earth’s  equator. 
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soid,  and  any  trajectory  path  that  came  lower  than  this  distance  was  considered  to  be  re-entrant  and  hence  for¬ 
bidden.  Geomagnetic  cutoff  rigidities  are  determined  by  calculating  charged  particle  trajectories  at  discrete 
rigidity  intervals  starting  with  a  rigidity  value  high  above  the  highest  possible  cutoff  and  decreasing  the  rigidity 
to  a  value  that  satisfied  our  criteria  that  the  lowest  allowed  trajectory  had  been  calculated.  As  the  calculations 
progress  down  through  the  rigidity  spectrum,  the  results  change  from  the  easily  allowed  orbits  to  a  complex 
structure  of  allowed,  forbidden,  and  quasi-trapped  orbits  (loosely  called  penumbra)  and  finally  to  a  set  of  rigidi¬ 
ties  where  all  trajectories  intersect  the  solid  earth.  Rigidity  interval  spacing  of  0.01  GV  was  used  throughout 
the  cosmic  ray  penumbra.  As  a  result  of  these  trajectory  calculations  we  determined  the  calculated  upper  cutoff 
rigidity  (Ru)  which  is  the  rigidity  value  of  the  highest  allowed/forbidden  pair  of  adjacent  cosmic  ray  trajecto¬ 
ries,  the  calculated  lowest  cutoff  rigidity  (Rl)  which  is  the  rigidity  value  of  the  lowest  allowed/forbidden  pair  of 
adjacent  cosmic  ray  trajectories,  and  an  ’’effective  cutoff  rigidity”  (Rc)  found  by  summing  the  allowed  orbits 
through  the  penumbra.  See  Cooke  et  al.  (1991)  for  definitions  of  cosmic  ray  cutoffs. 

Use  of  these  updated  three-hour  UT  interval  geomagnetic  cutoff  rigidity  models  results  in  better  fits  to  ex¬ 
perimental  data  than  the  previous  values  of  the  six-hour  UT  world  girds  (Smart  et  al.,  1999a,b,c).  A  compari¬ 
son  of  the  components  of  this  dynamic  cutoff  rigidity  model  with  vertical  cutoff  rigidities  calculated  using  only 
an  internal  geomagnetic  field  for  the  same  altitude  (Smart  and  Shea,  1997)  show  that  these  magnetospheric  cut¬ 
off  rigidity  values  are  consistently  lower  in  magnitude,  even  for  the  Kp=0  case. 

Rigidity  (momentum  per  unit  charge)  is  not  the  most  convenient  unit  for  use  with  energetic  particle  data 
since  most  energetic  particle  measurements  are  in  units  of  energy.  For  comparison  purposes,  we  have  selected 
the  invariant  latitude^  calculated  from  the  internal  geomagnetic  field  as  a  common  parameter.  We  have  interpo¬ 
lated  through  our  world  grids  of  vertical  geomagnetic  cutoff  rigidities  for  each  magnetic  activity  level  to  deter¬ 
mine  proton  cutoff  energy  contours  as  a  function  of  invariant  latitude  and  obtained  an  average  invariant  latitude 
for  each  energy.  The  invariant  latitudes  of  the  proton  cutoff  for  all  the  magnetic  activity  levels  are  illustrated  in 
Figure  1.  These  curves  indicate  an  almost  linear  relation  between  the  proton  cutoff  energy  with  magnetic  lati¬ 
tude  in  the  range  from  about  10  MeV  to  a  few  hundred  MeV.  We  note  that  the  change  of  proton  cutoff  energy 
with  Kp  is  relatively  uniform  over  the  range  of  the  original  Tsyganenko  (1989)  model  (Kp  values  0  through  5), 
but  the  cutoff  changes  introduced  by  the  Boberg  et  al.  (1995)  extension  are  non-linear  with  the  Dst  increment. 


PROTON  CUTOFF  ENERGY  (MeV) 

Fig.  1.  The  invariant  magnetic  latitude  of  the  proton 
effective  vertical  cutoff  energy  at  450  km  for  all  inte¬ 
ger  magnetic  activity  levels. 


Geomagnetic  cutoffs  are  angle-dependant.  Us¬ 
ing  Stormer  (1955)  theory  we  can  extrapolate  from 
the  vertical  cutoff  rigidity  to  the  cutoff  at  other  azi¬ 
muth  and  zenith  angles  by  the  following  equations: 

Rv  (1) 

Ra  =  4Rv  \  {[  1  +  (1  -  sin€  sincp  cos^A.)^']^  }.  (2) 

Rv  is  the  vertical  cutoff  rigidity,  Vk  is  the  constant, 
R«  is  the  cutoff  rigidity  in  a  specific  angular  direc¬ 
tion,  X  is  appropriate  magnetic  latitude,  e  is  the  an¬ 
gle  from  the  zenith  direction,  and  (p  is  the  azimuthal 
angle  measured  clockwise  from  magnetic  north. 
The  angular  geomagnetic  cutoffs  are  lowest  from 
the  magnetic  west  and  maximum  from  the  magnetic 
east.  In  Stormer  theory,  the  cutoffs  in  the  magnetic 
north  and  south  directions  are  the  same  as  the  verti¬ 
cal  cutoff 


*  Invariant  latitude  is  a  development  from  magnetospheric  physics.  It  is  a  magnetic  latitude  that  accounts  for  the  effective 
magnetic  center  displacement  from  the  geocenter,  the  orientation  of  the  magnetic  field  axis,  and  the  higher  order  terms  in 
the  magnetic  field. 
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ALTITUDE  INTERPOLATION 

The  basis  of  the  interpolation  technique  is  the  fact  that  the  Mcllwain  "L‘*  parameter  (Mcllwain,  1961)  can 
be  utilized  in  equation  (1)  for  cosine  squared  of  the  magnetic  latitude.  Then,  the  vertical  cutoff  rigidity  (Rc) 
equation  can  be  written  in  the  form 


Rc  =  Vk/Ll 


(3) 


The  latitude  and  longitude  of  the  location  for  which  we  want  to  determine  the  cutoff  rigidity  will  be  contained 
in  a  specific  "box”  on  the  5-degree  by  5-degree  world  grid.  Each  “box”  contains  the  Vr  and  L  parameters  for 
each  comer.  The  orientation  of  the  specified  latitude  and  longitude  position  within  the  grid  box  is  then  deter¬ 
mined.  We  use  the  L  value  of  the  latitude  and  longitude  (at  450  km)  and  the  Vr  constants  for  the  above  cutoff 
equation  at  each  "box  comer",  to  interpolate  a  cutoff  value  for  the  coordinate  (latitude  and  longitude)  in  the 
"box".  A  value  is  obtained  for  the  upper  cutoff,  Vr(ru),  the  lower  cutoff,  Vr(rl),  and  the  effective  cutoff,  Vr(rc). 
Then  the  Mcllwain  "L"  value  is  computed  at  the  specific  coordinate  and  for  both  altitudes  with  "L"  interpola¬ 
tion  used  to  extrapolate  from  the  values  at  450  km  altitude  to  the  specified  altitude. 


COMPARISON  WITH  SPACECRAFT  PROTON  CUTOFF  DATA 


In  order  to  test  the  accuracy  of  the  interpolated 
cutoff  values,  we  compare  our  computed  invariant 
latitude  of  the  solar  proton  cutoff  with  the  meas¬ 
ured  cutoff  latitudes  observed  by  the  SAMPEX 
spacecraft  (Leske  et  al.,1997).  The  solid  dark  line 
at  67.5  degrees  in  Figure  2  indicates  the  invariant 
cutoff  latitude  predicted  from  use  of  the  internal 
IGRF  magnetic  field  model  (Smart  and  Shea, 
1997).  The  gray  squares  indicate  the  computed 
cutoff  invariant  latitude  for  the  center  energy  of  the 
29-64  MeV  detector.  The  cutoff  latitude  for  each 
Kp  value  is  indicated  on  the  left  of  the  figure  and 
the  integer  Kp  values  are  on  the  right.  The  solid 
diamonds  are  the  cutoff  latitudes  of  Leske  et  al. 
(1997).  This  figure  shows  that  the  interpolated  and 
measured  cutoff  latitudes  exhibit  the  same  general 
trend,  for  both  magnetically  quiet  and  active  times. 


E  TSYGANENKO  MAGNETOSPHERIC  MODEL  PREDICTON  (Rc) 
♦  SAMPEX  29-64  MeV  PROTON  CUTOFF  MEASUREMENT 


DAY  OF  YEAR  1992 

Fig.  2.  Comparison  of  the  calculated  and  measured 
Invariant  latitudes  of  the  proton  cutoff  energy  as  a 
function  of  the  Kp  index. 


COMPARISON  WITH  SPACECRAFT  DOSIMETRY  DATA 

The  STS28  space  shuttle  flight  in  August  1989  was  in  a  57°  inclination  orbit  and  carried  a  dosimeter  that 
provided  data  on  the  large  solar  proton  event  sequence  that  began  on  12  August  19'89.  The  flight  was  com¬ 
pleted  on  13  August  with  the  shuttle  landing  shortly  after  13  UT.  We  used  the  dosimetry  data  acquired  during 
the  actual  space  flight  as  a  method  of  evaluating  the  accuracy  of  the  dynamic  cutoff  rigidity  model.  During  the 
time  the  vehicle  encountered  solar  protons,  the  geomagnetic  activity  level  as  quantified  by  the  Kp  index  varied 
between  1  and  3.  We  used  our  dynamic  cutoff  rigidity  model  to  predict  the  solar  proton  transmission  through 
the  magnetosphere  to  the  position  of  the  space  shuttle  (latitude,  longitude  and  altitude)  each  minute  of  the  shut¬ 
tle  fight  during  its  encounter  with  the  solar  proton  event.  The  solar  particle  flux  predicted  to  arrive  to  the  skin 
of  the  space  shuttle  was  used  as  input  to  the  NASA/JSC  PDOSE  code  to  compute  the  radiation  dose  rate.  A 
more  detailed  description  of  this  radiation  dose  rate  calculation  method  is  given  by  Golightly  and  Weyland 
(1997).  The  general  procedure  was: 
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1.  Determine  (by  interpolation  in  the  appropriate  5°  x  5°  world  grid  of  cutoff  rigidities)  the  vertical  cutoff  pro¬ 
ton  energy  at  the  spacecraft  position  for  each  minute  of  the  solar  proton  event  encounter. 

2.  Model  the  GOES  measured  solar  proton  flux  as  the  differential  flux  impacting  the  magnetosphere.  We  used 
a  three-section  fit  of  the  solar  proton  energy  spectrum:  an  exponential  between  30  and  50  MeV,  another  ex¬ 
ponential  between  50  and  60  MeV,  and  another  exponential  between  60  and  100  MeV.  The  fit  between  60 
and  100  MeV  was  used  to  extrapolate  to  higher  energies. 

3.  Determine  the  solar  proton  flux  at  energies  exceeding  the  proton  cutoff  energy  at  the  spacecraft  position. 

4.  Reduce  the  free  space  omni-directional  flux  to  account  for  “earth  shadowing”  in  low-earth  orbit. 

5.  Attenuate  the  incident  solar  proton  through  the  spacecraft  mass  distribution. 

6.  Calculate  the  dose  rate  from  the  attenuated  solar  proton  flux  spectrum  penetrating  to  the  dosimeter  location. 
Comparisons  Using  Cutoff  Parameters 

The  STS28  dose  rate  data  set  also  provided  an  opportunity  to  check  the  applicability  of  the  various  cutoff 
rigidity  parameters  such  as  Ru,  Rl,  and  Rc.  We  found  that  truncating  the  Kp  magnetic  activity  index  (i.e.  reduc¬ 
ing  the  +  sub  increment)  results  in  a  cutoff  value  at  the  spacecraft  position  that  is  too  high  and  the  solar  particle 
flux  predicted  to  be  allowed  at  the  spacecraft  was  not  sufficient  to  reproduce  the  dose  measurement.  It  is  neces¬ 
sary  to  “round  up”  the  Kp  index  in  order  to  obtain  a  cutoff  value  that  predicts  a  dose  rate  that  is  consistent  with 

the  measured  dose  rate.  , 

Evaluation  of  the  solar  proton  flux  allowed  to  the  spacecraft  using  each  of  the  of  three  cutoff  rigidity  pa¬ 
rameters  resulted  in  the  following: 

1 .  Use  of  the  upper  cutoff  rigidity,  Ru,  resulted  in  an  underprediction  of  the  dose  rate. 

2.  Use  of  the  lower  cutoff  rigidity,  Rl,  resulted  in  an  overprediction  of  the  dose  rate. 

3.  Use  of  the  effective  cutoff  rigidity,  Rc,  resulted  in  a  slight  systematic  underprediction  of  the  dose  rate. 

4.  Using  all  of  the  solar  particle  flux  down  to  Ru  and  then  using  the  average  transparency  of  the  penumbra  to 

specify  the  solar  particle  flux  transmitted  through  the  penumbra  between  Ru  and  Rl  resulted  in  a  systematically 
better  correspondence  between  the  computed  and  measured  dose  rate.  This  average  transparency,  a  function  of 
the  cutoff  rigidity,  is  shown  in  Figure  3. 


UPPER  CUTOFF  (Ru  in  GV) 
Fig.  3.  Average  transparency  of  the  penumbra. 


gust  1989  solar  proton  event. 


It  was  also  necessary  to  include  the  effect  of  the  lower  cutoffs  in  the  western  direction  to  reproduce  the  ob¬ 
served  radiation  dose  intensity/time  profile.  The  allowed  solid  angle  of  particle  access  was  divided  into  five 
equal  solid  angle  segments,  vertical,  north,  south,  east  and  west.  The  particle  flux  allowed  though  the  geomag¬ 
netic  cutoff  for  each  segment  was  summed  to  find  the  allowed  particle  flux  over  the  solid  angle  of  particle  ac¬ 
cess.  Employing  the  dynamic  cutoff  rigidity  model  with  a  proper  selection  of  the  magnetic  activity  index  re¬ 
sulted  in  a  minute-to-minute  correspondence  between  the  time  periods  of  computed  radiation  dose  rate  due  to 
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solar  protons  being  allowed  through  the  magnetosphere  to  the  position  of  the  space  shuttle  and  the  measured 
dose  rate  in  the  vehicle,  even  for  very  small  doses.  Employing  the  solar  proton  flux  spectral  and  the  appropriate 
Kp  magnetic  activity  indices  enabled  computation  of  the  solar  particle  dose  rate  for  each  minute  of  the  STS28 
flight.  The  comparison  between  the  computed  and  measured  dose  rate  is  given  in  Figure  4.  A  more  detailed 
description  of  these  results  are  given  by  Smart  and  Shea  (2002)  and  Smart  et  al.  (2003). 


EXTRAPOLATION  TO  THE  EARTH’S  ATMOSPHERE 


Since  we  have  a  method  of  extrapolating  verti¬ 
cal  cutoff  rigidities  to  various  spacecraft  altitudes,  it 
seems  apparent  that  we  can  use  the  same  procedure 
to  extrapolate  to  the  atmosphere.  We  will  demon¬ 
strate  this  by  calculating  vertical  cutoff  rigidities  for 
the  Concorde  flight  from  London  Heathrow  Airport 
(LHR)  to  New  York  Kennedy  airport  (JFK)  that  oc¬ 
curred  on  24  October  1989.  During  this  flight  high- 
energy  solar  cosmic  rays  (illustrated  in  Figure  5)  pro¬ 
vided  radiation  exposure  in  excess  of  the  background 
cosmic  radiation.  During  the  flight  time  the  Kp  mag¬ 
netic  index  was  4*^.  The  aircraft  position  coordinates 
provided  by  Dyer  (2002)  were  used  to  calculate  the 
vertical  cutoff  rigidity  and  the  vertical  proton  cutoff 
energy,  along  the  flight  path. 
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Fig.  5.  The  high-energy  solar  cosmic  ray  event  of  24 
October  1982  as  observed  by  high  latitude  neutron 
monitors.  Oulu,  Finland,  Goose  Bay,  Canada,  and 
Thule,  Greenland  are  positioned  near  the  boundaries 
of  the  North  Atlantic.  The  time  of  the  Concorde  flight 
is  indicated  by  the  box. 
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Since  this  was  an  interval  of  disturbed  geomag¬ 
netic  conditions  there  was  a  reduction  of  the  geo¬ 
magnetic  cutoff  below  the  normal  values.  Figure  6 
illustrates  the  normal  quiescent  cutoff  rigidity  values 
that  would  be  experienced  along  this  flight  path  to¬ 
gether  with  the  values  calculated  for  the  actual  geo¬ 
magnetic  conditions  for  this  specific  flight.  (Western 
hemisphere  longitudes  are  designated  as  negative  in 
this  Figure.)  Further  analysis  of  these  cutoff  rigidity 
data  and  the  corresponding  radiation  dose  experi¬ 
enced  at  flight  altitude  has  been  done  by  Dyer  et  al., 
(2003). 
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24  October  1989.  Bottom:  Calculated  vertical  proton 
cutoff  energy  along  the  fight  path  for  quiet  (Kp=0)  and 
actual  Kp=4  magnetic  activity. 
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FINAL  REMARKS 

We  have  developed  a  dynamic  geomagnetic  vertical  cutoff  rigidity  model  that  predicts  the  energetic 
charged  particle  transmission  through  the  magnetosphere.  This  reference  set  of  world  grids  covers  all  magnetic 
activity  levels  from  super  quiet  to  extremely  disturbed  (i.e.,  Kp  indiees  ranging  from  0  to  9+)  for  every  three 
hours  in  universal  time.  These  world  grids  of  vertical  cutoff  rigidities  were  obtained  by  particle  trajectory  trac¬ 
ing  in  model  magnetospheres.  Testing  of  the  uses  of  these  Interpolated  cutoff  rigidities  compares  very  favora¬ 
bly  with  actual  spacecraft  measurements  of  the  cutoff  latitude.  Use  of  these  cutoff  rigidity  models  for  comput¬ 
ing  radiation  dose  resulted  in  a  one-to-one  correspondence  between  the  portion  of  the  orbit  predicted  to  be  sub¬ 
jected  to  solar  protons  and  the  portion  of  the  orbit  where  the  radiation  dose  can  be  attributed  to  solar  particles. 
Although  developed  for  space  applications,  we  demonstrate  the  applicability  of  this  library  of  cutoff  rigidity 
models  for  computing  the  cutoff  rigidity  along  aircraft  routes. 
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